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ABSTRACT 

Heterogeneous clouds or temperature perturbations in rotating brown dwarfs produce variability 
in the observed flux. We report time-resolved simultaneous observations of the variable T6.5 brown 
dwarf 2MASSJ22282889-431026 over the wavelength ranges 1.1 - 1.7 /.tm and broadband 4.5 ^m. 
Spectroscopic observations were taken with Wide Field Camera 3 on board the Hubble Space Telescope 
and photometry with the Spitzer Space Telescope. The object shows sinusoidal infrared variability 
with a period of 1.4 hr at most wavelengths with peak-to-peak amplitudes between 1.45% and 5.3% 
of the mean flux. While the light curve shapes are similar at all wavelengths, their phases differ 
from wavelength to wavelength with a maximum difference of more than half of a rotational period. 
We compare the spectra with atmospheric models of different cloud prescriptions, from which we 
determine the pressure levels probed at different wavelengths. We find that the phase lag increases 
with decreasing pressure level, or higher altitude. We discuss a number of plausible scenarios that 
could cause this trend of light curve phase with probed pressure level. These observations are the first 
to probe heterogeneity in an ultracool atmosphere in both horizontal and vertical directions, and thus 
are an ideal test case for realistic three dimensional simulations of the atmospheric structure with 
clouds in brown dwarfs and extrasolar planets. 

Subject headings: brown dwarfs — stars: atmospheres — stars: individual 
(2MASSJ22282889-4310262) — stars: variables: general 



1. INTRODUCTION 

Flux variability in brown dwarfs has been attributed to 
heterogeneous cloud cover that modulates the light curve 
through rotation. Two early T dwarfs have recently been 
discovered as photometrically va r iable at the multipl e 
percent level ()Artigau et all 120091: iRadigan et aLll2012[ ). 
These objects lie at the L/T transition between cloudy L 
dwarfs and mostly clear late T dwarfs, where clouds have 
dissipated or sunk below the photosphere. The mecha- 
nism for cloud dispersion is not well understood, and dif- 
ferent models have been proposed that invoke cloud thin- 
ning and rain-out (e g.. iKnapp et al.|[200l iTsuiil 120051: 
iBurrows et all 120061: ISaumon fc Marlevll2008D or patchv 
clouds (jMarlev et al.l I2010| r Variability indicates het- 
erogeneous cloud structure, but models with cloud holes 
that allow flux from hot deep regions to escape do not 
fit the measured time-resolved spectroscopy in early T- 
dwarfs (jApai et al.l Isubmittedl ) . A mixture of thin and 
thick clouds provides a better match. 

Variability is not limited to L/T transition dwarfs, 
though these appear to show the highest variability 
level. I band surveys found seve r al variable L dwarf s 
(e.g. iBailer- Jon es fc Muiidtl [200lt iGelino et al.] [200l . 
and iClarke et al. (2008) reported variability for two T 
dwarfs beyond T5. While later T dwarf spectra have 
been modeled as clear atmospheres with reasonable suc- 
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cess (e.g. lAllard et al.l 120031: IStephens et al.ll2009[ ). color 
dispersion and variability indicat e that clouds, such as 
sulfide clouds (jMorlev et al.l 120121 ) . may persist in these 
atmospheres. 

The T6.5 dwarf 2MASSJ22282889-431 0262 (h ereafter 
2M2228) was disc overed by Burgasser 'et al.l ([2003). 
IClarke et al.l (|2008D measured periodic variability, likely 
due to rotation with a period of 1.43 hr and a peak- 
to- valley amplitude of 1.4% ± 0.14% of the mean flux 
in J band. Its pa.rallax was measured as 94 ± 7 mas 
by iFahertv et all H2012D . resuhing in a distance d — 
10.64 ±0.79 pc. 

Here we report time resolved high-precision observa- 
tions of 2M2228: simultaneous HST/^NYCi spectroscopy 
and Spitzer/lYlKC photometry that for the flrst time si- 
multaneously probe variability at different altitudes in 
the atmosphere of a brown dwarf. The observations re- 
veal a phase shift between wavelengths that is evidence 
for vertical heterogeneity in the atmosphere. They pro- 
vide unique constraints on the heterogeneous cloud or 
temperature structure in both the horizontal and verti- 
cal directions. 

2. OBSERVATIONS 

We observed 2M2228 on 2011 July 07 with the in- 
frared channel of the Wide Field Camera 3 (WFC3) on 
the Hubble Space Telescope (HST) and the Infrared Ar- 
ray Camera (IRAC) of the Spitzer Space Telescope (GO- 
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12314, GO-70170; PI: D. Apai). T/ST observed between 
BMJD(UTC) 55749.528 to 55749.886, Spitzer between 
55749.766 and 55750.002, providing an overlap of almost 
three hours. 

With HST, we obtained spectral time series with the 
G141 grism for the wavelength range 1.05 — 1.7 /zm over 
six consecutive HST orbits, spanning 8.58 hr. For each 
orbit of 96 minutes, we observed the target for 35.6 min- 
utes, obtaining nine exposures of 224 s. For the rest of 
the time, the target was occulted by Earth. 

The infrared channel of WFC3 is a Teledyne HgCdTe 
detector with 1024x1024 pixels. We used the 256 x 256 
pixels subarray with a field of view of approximately 30 
X 30"to avoid WFC3 buffer dump overheads. Each expo- 
sure is constructed from multiple non-destructive read- 
outs. We used the SPARS25 readout mode with 11 reads 
per exposure. Each sub-readout has an exposure time of 
22.34 s after an initial first short read at 0.27 s. The 
maximum number of counts in a pixel at the end of an 
exposure was 27,000 counts, well below the detector half- 
well capacity. 

The first-order spectrum has a width of ~ 140 pixels 
and a dispersion of 4.65 nm pixel"^. The zeroth and 
second orders were not recorded. A direct image was 
taken for wavelength calibration through the F132N nar- 
row band filter at the beginning of each orbit. We kept 
the spectra on the same pixels for all orbits to avoid sys- 
tematic errors from imperfectly corrected pixel-to-pixel 
sensitivity variations. 

The Spitzer observations were taken with IRAC chan- 
nel 2 at 4.5 ^m (warm mission phase). Channel 2 has a 
256 X 256 pixel InSb detector with a field of view of 5.2 x 
5.2 arcmin. We acquired 2876 individual measurements 
with 4.4 s exposure time and cadence of 7.2 s, covering 
5.75 hr. 

3. DATA REDUCTION 
3.1. HST/WFC3 

We based our data reduction on the combination of 
standard WFC3 pipeline, the PyRAF software package 
aXe0, and custom IDL routines. We summarize here the 
m ost important poin ts, additional details can be found 
in lApai et al.l ([submitted ). 

The WFC3 pipeline calf w3 produces flux calibrated, 
two-dimensional spectral images with flagged bad pixels. 
We start our reduction from the . ima files which contain 
an array of each subread for each exposure. We discard 
the short zero read, resulting in 10 images of 22.34 s per 
exposure. After the full reduction and spectral extrac- 
tion, we average the 10 spectra obtained per exposure 
into one to reduce the noise. 

We correct bad pixels by interpolating over neighbors 
and cosmic rays by replacing outliers in a subexposure 
sequence with the median pixel value over the sequence. 
For the background subtraction, flat fielding and extrac- 
tion of the two dimensional spectra, we use aXe,^ de- 
veloped specifically for HST slitless spectroscopy. We 
choose the spectral extraction width to be three times 
the FWHM of the vertical spectral profile. The FWHM 
is determined individually for each orbit by averaging all 
spectra and collapsing the result into a one-dimensional 
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(ID) vertical profile, to which we fit a Gaussian. The 
mean extraction width was 5.651 ±0.013 pixels. The po- 
sitional differences of the spectra over a visit were deter- 
mined with cross-correlation to be smaller than 0.02 pix- 
els. 

The extracted ID spectra are flux calibrated with the 
G141 sensitivit y curve, which has an estimated uncer- 
tainty of 1% (|Kuntschner et al.l 1201 ID . We calculate 
the flux uncertainty from the read-out noise and the 
photon noise of the source and background. The time 
of observation is calculated as the mid-time between 
sub-reads from the header keywords "ROUTTIME" and 
"EXPTIME". The times are converted from Modified 
Julian Date (MJD) to Barycentric Modified Julian Date 
(BMJD-UTC) with the barycen.pro IDL routine to allow 
comparison with the Spitzer data, and double-checked 
with the utc2bjd.pro routine. These routines are accu- 
rate to well below a second (|Eastman et al.ii201Qi . The 
spectra have a useable wavelength range of 1.05 — 1.7 /^m 
and a resolution R = 130. 

We identified and corrected two detector systematics 
at sub-percent levels. First, in each batch of sub-reads 
of an exposure the measured flux dropped by 1% — 2% 
from the first to last su-bread for the bright pixels in a 
spectrum, likely a result of charge diffusion. A second 
effect was an exponential ramp at the beginning of each 
orbit, which was found to be independent of count rate 
and wavelength. Both effects were co rrected with an- 
alytic functions (jApai et al.l Isubmittedl ) derived from a 
non- variable source (ramp) or an average over all sources 
(fiux drop). 

3.2. Spitzer /IRAC 

The raw Spitzer data were processed with the IRAC 
pipeline ver. 19.0.0, which performs dark subtraction, 
fiat-fielding, fiagging of bad pixels and flux calibration. 
We masked columns that were significantly noisier than 
the rest of the pixels, and applied a position-dependent 
fiux correction. We determined the source center with 
the gcntr.pro routine with an FWHM of 4 pixels. We 
performed aperture photometry with the aper.pro rou- 
tine with an aperture of 2 pixels radius, and the sky 
annulus between 12 and 20 pixels radius, which mini- 
mized the noise. We identified and removed frames with 
cosmic ray hits across the source by checking for anoma- 
lous centroid positions and flux outliers of more than 5cr. 
To correct for intrapixe l sensitivity variations, we used a 
pixel mapp ing method (jKnutson et al.l[20T^ iLewis et"all 
Isubmittedl ). We determined the noise pixel parameter 

/3 for each exposure and created a pixel flux map using 
each measurement's 50 nearest spatial neighbors. After 
applying the flux correction, we binned the data into 5 
minute bins to reduce the noise. We calculated the er- 
ror from the shot noise and the statistical uncertainty of 
the sky subtraction. The remaining systematic errors af- 
ter the intrapixel sensitivity correction were smaller than 
the random errors in our binned data. The time of ob- 
servations in BMJD(UTC) is given in the Spitzer data 
header BMJD_OBS, and was adjusted to correspond to 
the mid-exposure time. 

4. RESULTS 

2M2228 shows periodic variability in both the HST 
and Spitzer daia, sets. Figure [1] (upper panels) shows the 
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Fig. 1. — HST (blue) and Spitzer (red) observations of 2M2228 integrated over different bandpasses. Top: ffST time series integrated 
over 1.1 — 1.7 /im and Spitzer pfiotometry at 4.5 /im as a function of barycentric time. Bottom: period folded ligfit curves in six wavelength 
regions folded with the rotation period P = 1.405 hr. A representative error bar at each wavelength is a combination of random and 
systematic errors. The black curve is the best sine fit. The phase is given with respect to the phase of the narrow J light curve. 

HST time series integrated over the full wavelength range 
compared to the Spitzer photometry. 2M2228 is a fast 
rotator with a rotation period of only 1.405 ± 0.005 hr 
determined from a sine fit to the HST data. Each HST 
orbit of 96 minute samples part of a subsequent rota- 
tion period with a shift by l/6th of a period. Six orbits 
cover the full rotation phase, and 90% of the phases are 
sampled at least twice. The Spitzer observations contin- 
uously cover four adjacent rotation periods, two of which 
partially overlap in time with the HST data. The derived 
period is 1.38 ±0.03 hr. We adopt the more precise HST 
period for further analysis. 

The mean HST spectrum is plotted in Figure [2] (top 
panel) together with a ground-based SPEX spectrum 
(|Burgasser et al.|[2004D and compared to models (see Sec- 
tion [S|) . The spectra are dominated by deep water ab- 
sorption features with higher flux in the J and H band 
windows. Other features include CH4 and K I and Na 



I. We study the variability as a function of wavelength 
by integrating the spectra for a number of characteristic 
regions. The bandpasses are selected such that the flux 
is emerging from a well-deflned pressure region in the at- 
mosphere. Integrated and period-folded light curves at 
these wavelengths, together with the Spitzer photometry 
at 4.5 /im, are shown in the lower panels of Figure [T] 

The light curve shapes in all wavelength regions are 
well-flt by sinusoidal functions of the form F{t) — 1 + 
Aq cos t + Bosint. We determine amplitudes, phases and 
errors from sine flts by bootstrapping (Figures [1] and [J) . 
We do not find significant variability for the wavelength 
range 1.12 — 1.17. Because of high noise due to low grism 
sensitivity and poor correction of systematics there may 
be undetected variability on the l%-2% level. We do 
not include this wavelength range in our further analy- 
sis. The peak-to-valley amplitude levels range between 
« 1.5% in the IRAC band more than 5% in the deep water 
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band around 1.4 /im. 

Most interestingly, we find significant shifts in the 
phase of the hght curve as a function of wavelength. This 
is evident in the lower panels of Figure [1] All light curves 
are period- folded with reference to the same time tg , tak- 
ing the narrow J band light curve as zero reference for 
the phase. The narrow H band curve shows a small delay 
in phase of 14°.6 ± 2°, or 3.4 ± 0.5 minutes. The 4.5 
phase curve is delayed much more: 118±7° or 27.7±1.6 
minutes. The deep absorption bands show the largest 
phase difference to J: 165° ± 7° for 1.62 — 1.69 /im and 
189° ± 6° for 1.35 — 1.43 /im, more than half a rotation 
period. 

5. MODELING 

Although current state-of-the-art atmospheric models 
of brown dwarfs are limited to a ID treatment of ra- 
diative transfer and convection, and do not include local 
heterogeneities, they can provide globally averaged atmo- 
spheric quantities. We fit ID models to the mean spec- 
trum, and to the SPEX spectrum from iBurgasser et al.l 
(|2004f ) for broader coverage into the Y and K bands, 
to derive basic atmospheric parameters and determine 
the pressure levels that are probed by different wave- 
lengths. We use two independent sets of T-dwarf mod- 
els to test the robustness of our conclusions. The 



first are models by 


Morlev 


et al.l 


lAckerman & Marlev 


|200l|) 


cloud 



sulfide and chromium clouds. The models use a sedi- 
mentation efficiency parameter fsed, where larger fsed 
corresponds to larg er grains and th i nner c louds. The sec- 



ond are mode ls bv iBurrows et al. (|20Q6|) a nd adaptions 
thereo f froi n iHubenv &: Burrows (I2007D . iCurrie et al.l 
(poll . and lMadhusudhan et all (|2011h that include dif 



ferent chemical compositions, cloud thickness prescrip- 
tions and grain sizes. 

We find the best spectral match using a Morley model 
(Figure[21) with effective temperature T^ff = 900 K , sur- 
face gravity of log 5 = 5, and fsed = 5, though \ogg — 4.5 
and fsed = 3 (thicker clouds) provide an almost equally 
good fit. This model includes thin Na2S, MnS, and Cr 
clouds. The overall spectral fit is good, with the main 
discrepancy in the H band around 1. 65 nm.. This discrep - 
ancy is common m many sources (IMorlev et al.l[20Tl . 
and is likely the result of incomplete methane opacities. 
A second discrepancy are alkali absorption features at 
1.24 /im and 1.12 /im, whose depths are overestimated 
by all models. 

The best Burrows model is an "iron" cloud model with 
Teff = 900 K and logg = 4.5, where iron provides rep- 
resentative optical properties. It uses the "Model A" 
cloud prescription (Burrows ct al. 200f|), where the cloud 
density is exponentially declining towards the top. The 
overall spectral shape is reproduced reasonably well, al- 
though with similar and somewhat larger discrepancies 
than the Morley model. Again, other models with differ- 
ent cloud prescriptions or chemical abundances provide 
similar fits. We find that the main conclusions of this 
work do not depend on the detailed cloud model struc- 
ture. 

From our best-fit models we infer the pressure levels 
that are typically probed at different wavelengths. With 
this information we translate the phase shift in the light 
curves as a function of wavelength to a phase shift as 



a function of pressure. Together with the longitudinal 
information from the light curve variability, our observa- 
tions constrain heterogeneity both in the horizontal and 
vertical direction. 

Figure [2] (bottom) shows the pressure levels at optical 
depth T — 2/3, with cloud levels for the Morley model in- 
dicated. In reality, a range of pressures is probed simulta- 
neously, a nd can be described by contribution functions 
(see, e.g., iMarlev et al.]|2012|) . Instead of using these 
functions, because we are only looking for trends as a 
function of probed pressure region, for each wavelength 
band we derive an average pressure level and assign half 
an atmospheric scale height for the width of the region. 

In Figure [3] (left), the phase shift is shown as a func- 
tion of pressure for both models. We find that the phase 
delay increases strongly with decreasing pressure (higher 
altitude) for the four bandpasses where the model fit is 
good. For the 1.62-1.69 /^m bandpass, the model fit is 
poor because the methane opacities are likely underes- 
timated. This means that the actual probed pressure is 
lower (higher altitude) than determined from the model, 
making that point consistent with the trend. Qualita- 
tively, the trend is the same for both models, only shifted 
to lower pressures for the Burrows model, mostly because 
of the lower gravity of the model. 

The right panel of Figure [3] gives the peak-to-peak am- 
plitude of the variability as a function of pressure. There 
is no obvious correlation between amplitude and probed 
pressure level. 

6. DISCUSSION 

We have discovered sinusoidal infrared variability in a 
late T-dwarf with phase offsets between different wave- 
lengths. From comparing atmosphere models to the spec- 
trum, we find a clear trend of the phase shift with the 
pressure probed at respective wavelengths: the lower the 
pressure of the atmospheric layer, the larger the phase 
lag with respect to the highest pressure layer, up to half 
a rotation period. We consider three cases that may 
cause infrared variability to be out of phase: (I) change 
in the temperature-pressure profile T(p) without change 
in opacities, (II) change in the opacities (clouds and/or 
gas) without change in T{p), or (III) an intermediate 
scenario. 

Case (I): Define AT as the horizontal temperature dif- 
ference on isobars. To get the observed variability re- 
quires that regions where AT is positive aloft exhibit neg- 
ative AT at depth, and vice versa. Hotter regions may 
result from subsidence, which advects high-entropy gas 
from above, while cooler regions may result from ascent, 
which advects low-entropy gas from b elow. This s uggests 
a "stacked cell" scenario (e.g. Fletcher et al.l[20TTI ) where, 
on one side of the brown dwarf, subsidence at depth ac- 
companies ascending motion aloft — and vice versa on the 
other side. 

Case (II): In this scenario, changes to the brightness 
temperature result from the r = 2/3 level moving up 
or down. It seems unlikely that opacity changes due 
to clouds alone can produce the observed trend for the 
phases of the light curves. The strongest variability, 5%, 
is seen at the lowest pressure level, where no significant 
cloud opacity is expected from the models. If present, 
cloud opacity variations at such high altitude would af- 
fect the emission in other bands and cause the variability 
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Fig . 2. — Top: normalized mean HST/WFCS spectrum of 2M2228 together with the SPEX spectrum of the source by IBurgasser et aLl 
II2004I ). and best-fit atmospheric models by Morley and Burrows (see text). Bottom: pressure level probed at optical depth r = 2/3 as a 
function of wavelength for the two models. Cloud layers are indicated for the Morley model. Selected bandpasses to study the variability 
are indicated by the horizontal green lines (cf. Figures [l]and|3ll, and averaged pressures in these bandpasses are marked with dots. 
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Fig. 3. — Phase shifts (left) and amplitudes (right) of the observed 2M2228 light curves as a function of the dominantly probed pressure. 
There is a correlation between phase shift and pressure (shaded regions), but not between amplitude and pressure. Phase shifts in degrees 
are relative to the J band light curve. The pressure in each wavelength region is derived from the two models presented in Figure [2] The 
assumed width of the pressure region is half a scale height, the horizontal width of the bar is the ilc observational error. The pressure for 
bandpass 4 is likely overestimated due to underestimated methane opacity in the model and should be at lower pressure (arrow), making 
it consistent with the trend of the phase shift. 
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to be mostly in phase. A second, lower level cloud deck 
would need to show much stronger spatial variations in 
opacity of opposite phase to compensate for the opacity 
variations of the high cloud layer. 

Gas opacity is strongly wavelength dependent. In prin- 
ciple, regions with increased gas opacity at 1.4 /im could 
have decreased opacity in J and H, resulting in a phase 
difference in the light curves between these wavelengths. 
Still, pure gas opacity changes are unlikely since clouds 
likely still play a role in these atmospheres as found 
from the model comparison. A mixed scenario could 
be plausible, where the variability at higher pressure 
regions stems from changes in cloud opacity, while at 
low pressure the gas opacity varies with different phase. 
It is not clear which gas could produce such opacity 
changes. Water does not condense at the given tem- 
peratures, and C O and CII4 should largely be quenched 
(jHubeny h Burrows 2007) , making spatial variations dif- 
ficult. A condensin i g species, e.g. those considered in the 
iMorlev et al.l ()2012l ) models, where condensation induces 
spatial variations in the gas abundance, is a possibility. 

Case (III): Of various possibilities, a plausible scenario 
involves a heterogeneous deep cloud controlling the vari- 
ability at J and H, with horizontal variations in tem- 
perature aloft producing the variations in bands sensing 
high altitudes. Again, a "stacked cell" scenario as in case 
(I) would be a possible cause, where a region of subsi- 
dence that affects the deep cloud is overlain by a region of 
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ascent that affects the high-altitude temperature. Such 
a stacked-cell scenari o has been suggested for the tro- 
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